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THE  EFFECT  OF  CONSTITUTIVE  LAW  PERTURBATIONS 
ON  FINITE  ANTIPLANE  SHEAR  DEFORMATIONS 
OF  A  SEMI-INFINITE  STRIP 

By 

C.  O.  HORGAN  {University  of  Virginia,  Charlottesville,  Virginia) 

AND 

L.  E.  PAYNE  (Cornell  University,  Ithaca,  New  York ) 

Abstract.  This  paper  is  concerned  with  assessing  the  effects  of  small  perturbations 
in  the  constitutive  laws  on  antiplane  shear  deformation  fields  arising  in  the  theory 
of  nonlinear  elasticity.  The  mathematical  problem  is  governed  by  a  second-order 
quasilinear  partial  differential  equation  in  divergence  form.  Dirichlet  (or  Neumann) 
boundary-value  problems  on  a  semi-infinite  strip,  with  nonzero  data  on  one  end  only, 
are  considered.  Such  problems  arise  in  investigation  of  Saint- Venant  end  effects  in 
elasticity  theory.  The  main  result  provides  a  comparison  between  two  solutions,  one 
of  which  is  a  solution  to  a  simpler  equation,  for  example  Laplace’s  equation.  Three 
examples  involving  perturbations  of  power-law  material  models  are  used  to  illustrate 
the  results. 

1.  Introduction.  The  equilibrium  equations  governing  finite  antiplane  shear  defor¬ 
mations  of  some  homogeneous  isotropic  compressible  or  incompressible  nonlinearly 
elastic  materials  have  been  shown  to  reduce  to  a  single  second-order  quasilinear  par¬ 
tial  differential  equation  in  two  independent  variables  for  the  out-of-plane  displace¬ 
ment  (see  e.g.  [1-4]).  In  particular,  for  the  generalized  neo-Hookean  incompressible 
materials  for  which  the  strain-energy  depends  only  on  the  first  invariant,  the  gov¬ 
erning  equation,  in  the  absence  of  body  forces  is  (1.1)  below.  This  equation  also 
governs  finite  antiplane  shear  deformations  for  a  certain  class  of  compressible  mate¬ 
rials  [4j.  Considerable  attention  has  been  paid  to  the  analysis  of  solutions  of  (1.1) 
on  rectangular  domains  whose  lengths  greatly  exceed  their  widths.  In  particular,  for 
such  long  thin  domains  (or  for  semi-infinite  strips),  with  traction-free  lateral  sides, 
the  asymptotic  behavior  of  solutions  of  (1.1),  as  the  axial  variable  increases,  is  of 
interest  in  connection  with  Saint-Venant’s  principle  (see  e.g.  [5-8]  and  the  references 
cited  therein.  Recent  reviews  on  Saint-Venant’s  principle  are  given  in  [9,  10].  See 
also  [11]  for  consideration  of  a  fourth-order  analog  of  (1.1).)  The  spatial  evolution 
of  solutions  of  (l.l)  for  both  Dirichlet  and  Neumann  boundary  conditions  on  the 
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lateral  sides  is  also  relevant  to  the  study  of  Phragmen-Lindelof  type  principles  (for 
such  results,  in  two  or  three  dimensions,  see  e.g.  [12,  13]  and  the  references  cited 
therein).  Recently,  the  present  authors  have  investigated  the  asymptotic  behavior  of 
solutions  of  inhomogeneous  equations  which  are  generalizations  of  (1.1).  In  [14],  in¬ 
homogeneous  equations  where  a  constant  term  is  added  to  the  left-hand  side  of  ( 1 . 1 ) 
have  been  studied.  In  the  antiplane  shear  context,  this  would  correspond  to  a  constant 
body-force.  It  was  shown  in  [14]  that  solutions  to  Dirichlet  problems  for  such  equa¬ 
tions  are  well  approximated,  away  from  the  ends  of  a  finite  rectangle,  by  solutions 
to  the  corresponding  one-dimensional  problem  for  an  ordinary  differential  equation 
on  the  cross  section  of  the  rectangle.  Such  results  are  of  interest  in  assessing  the 
approximate  nature  of  one-dimensional  theories  compared  to  exact  two-dimensional 
theories,  and  have  played  an  important  role,  for  example,  in  establishing  plate  and 
shell  theories  in  solid  mechanics.  Applications  to  problems  in  geometry,  for  example, 
to  the  equation  of  a  surface  of  constant  mean  curvature,  have  also  been  discussed  in 
[14].  Generalizations  of  these  results  to  the  equations  governing  capillary  surfaces 
and  extensible  films  are  described  in  [IS], 

In  this  paper,  we  return  to  the  homogeneous  equation  (1.1)  on  a  semi-infinite 
strip,  subject  to  nonzero  Dirichlet  or  Neumann  boundary  conditions  on  the  near  end 
only,  and  examine  the  effect  of  perturbations  of  the  coefficient  p  on  solutions.  In 
the  context  of  antiplane  shear  [1-4],  this  coefficient  is  the  derivative  of  the  strain- 
energy  density  and  so  our  concern  is  with  the  effect  of  constitutive  law  perturbations 
on  solutions.  Such  results  are  of  interest  given  the  practical  difficulty  in  constructing 
constitutive  models  that  provide  an  exact  description  of  material  behavior. 

Specifically,  we  are  concerned  with  second-order  quasilinear  partial  differential 
equations  in  two  independent  variables  of  the  form 

=  0,  q  =  («  fiu  p)1'2 ,  (1.1) 

where  the  usual  summation  convention  is  employed  with  subscripts  preceded  by  a 
comma  denoting  partial  differentiation  with  respect  to  the  corresponding  Cartesian 
coordinate.  As  mentioned  above,  in  the  context  of  antiplane  shear,  p  is  determined 
by  the  constitutive  model  governing  material  behavior;  and  in  (1.1)  u  =  u(xl ,  x2)  is 
the  displacement  field.  A  commonly  used  constitutive  law  gives  rise  to  functions  p 
of  power-law  form, 

p  =  p(l  +bq2/n)"~' ,  p,b,n>  0,  (1.2) 

where  p  is  the  shear  modulus  for  infinitesimal  deformations,  and  n  is  a  material 
hardening  parameter.  The  case  n  =  1  in  (1.2)  corresponds  to  the  neo-Hookcan 
material  for  which  p  is  a  constant,  and  (1.1)  is  Laplace’s  equation.  Letting  p  = 
dp/dq2 ,  we  see  from  (1.2)  that  p  >  0  or  p  <  0  according  as  n  >  1  or  n  <  1 
respectively.  The  material  is  said  to  harden  or  soften  in  shear  in  these  situations 
[5].  When  n  =  1/2  in  (1.2),  equation  (1.1)  is  reducible,  by  a  change  of  scale,  to  the 
minimal  surface  equation 

(1  +«22)«  ,,  -  2«  ,«  2U  ,12  +  (1  +M2|)M.22  =  °- 


(1.3) 
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We  consider  both  Dirichlet  and  Neumann  problems  for  ( 1 . 1 )  on  the  semi-infinite 
strip  R  =  {(jc, ,  jc2)|0  <  x2  <  h ,  jc,  >  0} .  For  ease  of  exposition,  we  confine 
attention  in  what  follows  to  the  Dirichlet  problem.  The  modifications  necessary  to 
treat  the  corresponding  Neumann  problem  are  described  in  Sec.  4.  We  assume  the 
existence  of  classical  solutions  u  e  C2{R)  n  C'(R)  satisfying  (1.1)  on  R  subject  to 
the  boundary  conditions 

u{xt ,  0)  =  0,  u(xt,h)  =  0,  x,  >0,  (1.4) 

u,  u  ,  -*  0  (uniformly  in  x2)  as  jc, -♦  oo,  (1.5) 

u(0,  x2)  =  /(x2) ,  0  <  x2<  h,  (1.6) 

where  the  prescribed  function  /  is  sufficiently  smooth  and  satisfies  /(0)  =  f(h)  =  0 . 

For  a  rather  general  class  of  functions  p  it  was  established  in  [7]  that  solutions 
of  (1.1),  (1.4)— (1.6)  decay  exponentially  with  the  distance  jc,  from  the  end  x,  =  0. 
The  exponential  decay  rate  was  characterized  explicitly  in  [7]  in  terms  of  the  function 
p  and  the  strip  width  h  .  The  hypotheses  made  in  [7]  concerning  p ,  which  we  shall 
also  assume  here,  are  conveniently  separated  into  two  cases.  It  is  assumed  that  there 
exist  positive  constants  ma ,  Ma  and  nonnegative  constants  Ka  (a  =  1 , 2)  such 
that,  for  all  solutions  u  of  (1.1),  (1.4)— (1.6)  on  R,  we  have  either 
Case  1. 

0  <  m,  <  p  <  Af,  +  K,q2p,  (1.7) 

or 

Case  2. 

0  <  m2<  p~l  <  M2  + K2q2p,  (1.8) 

respectively.  As  pointed  out  in  [7],  if  p  were  a  bounded  function  of  its  arguments, 
then  the  AT,  in  (1.7)  could  be  taken  to  be  zero.  Roughly  speaking,  the  first  term  on 
the  right  in  (1.7)  provides  a  bound  on  p  as  q  — ♦  0,  while  the  second  term  gives  a 
bounding  function  for  p  as  q  oo .  A  function  p ,  for  which  (1.7)  holds,  is 

p  =  p(l+q2),  (1.9) 

which  may  be  viewed  as  a  special  case  of  (1.2)  with  n  =  2,  b  =  2.  For  this  p ,  we 
can  take  ml  =  p,  Mx  =  p ,  Kx  -  1  in  (1.7).  For 

p  =  p(l  +  2bq2)~1'2 ,  (1.10) 

(corresponding  to  the  value  n  =  1/2  in  ( 1 .2)),  in  which  case  (1.1)  is  reducible  to  the 
minimal  surface  equation  (1.3),  then  (1.8)  is  satisfied  with  m2  =  p~l ,  M2  =  , 

and  K2  =  2 bp~2 .  We  observe,  as  in  [7],  that  neither  (1.7)  nor  (1.8)  requires  that 
equation  (1.1)  be  elliptic,  that  is,  for  all  solutions  u  and  at  all  points  of  R , 

p  +  2p  q2  >  0  (p'  =  dp/dq2),  (1.11) 

although  the  results  obtained  in  [7]  and  here  are  primarily  of  interest  for  elliptic 
equations.  It  should  also  be  noted  that,  in  choosing  the  constants  mn  ,  Mn  in  ( 1 .7), 
(1.8)  for  a  given  p ,  it  is  desirable  to  choose  mn  as  large  as  possible  and  Ma  as 
small  as  possible. 
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In  this  paper,  we  are  concerned  with  comparing  solutions  u  of  (1.1),  ( 1 .4)— ( 1 .6) 
to  solutions  v  of  a  differential  equation  arising  from  a  small  perturbation  of  the 
constitutive  function  p .  The  comparison  solution  v  satisfies  the  same  boundary 
conditions  as  u .  In  particular,  such  a  comparison  is  of  interest  where  the  problem 
for  t;  is  much  simpler  than  that  for  u ,  for  example,  v  may  be  the  solution  to  a 
linear  problem.  Indeed,  for  two  of  the  illustrative  examples  considered  in  this  paper, 
v  is  a  harmonic  function.  We  assume  the  existence  of  v  e  C2(R)nCl  (R)  satisfying 


[P(p2)vj.a  =  0,  P  =  (v,fiV'fi)'/2,  (1.12) 

on  R ,  subject  to  the  boundary  conditions 

t>(x,,0)  =  0,  v(xl,h)  =  0,  *,>0,  (1.13) 

v,  v  ,  ->  0  (uniformly  in  x2)  as  x,  -*  oo,  (1.14) 

v(0,x2)  =  f(x2),  0  <  x2<  h.  (1.15) 

The  perturbed  constitutive  function  p  is  such  that  it  satisfies  conditions  analogous 
to  (1.7),  (1.8),  namely, 

0<ml<p<Ml+Klp2p  (1.16) 

or  _ 

0<m2<{p)~l  <M2  +  K2p2p.  (1.17) 


Thus,  in  Case  1,  both  (1.7)  and  (1.16)  are  assumed  to  hold  while  in  Case  2,  both 
(1.8)  and  (1.17)  are  to  hold.  (As  was  remarked  in  connection  with  (1.7),  (1.8),  it 
is  desirable  to  choose  ma  as  large  as  possible  and  Ma  as  small  as  possible.)  For 
e  «  1 ,  it  is  also  assumed  that  there  exist  positive  functions  y,(p ,  e) ,  y2(p ,  e) ,  with 

yx(j>,  e)  <c,(e)  <  1,  (1.18) 

V2(P,e)<c2 ,  (U9) 

where  c2  is  a  constant  such  that 

P\p(d2)  ~  P(P2)\  <y,iP,  e)p(q2)\q  -p\  +  e2 py2{p ,  e)[p(q2)p(p2)]'r~ .  (1.20) 

The  hypothesis  (1.20),  though  somewhat  complicated,  serves  to  define  the  consti¬ 
tutive  law  perturbation  of  concern  here.  The  form  of  (1.20)  arises  from  the  following 
considerations.  We  write 

p(d2)-p(p2)  =  (l/2)[p(q2)-p(p2)  +  p(p2)-p(p2)]  2 

+  (l/2)[p(q2)  -  p(q2)  +  p(q2)  -  p(p2)) 

so  that 

I P(q2)  ~  P(P2) I  <  (1/2)[| p(q2)  -  p(p2)\  +  | p(q2)  -  P(P2) II  22 

+  (1/2)[| p(q2)  -  p(q2) |  4- 1 p{p2)  -  p(p2) |] . 

The  second  two  terms  on  the  right  in  (1.22)  measure  the  difference  between  p(s) 
and  p(s)  at  the  same  value  of  their  arguments  while  the  first  pair  of  terms  measure 
the  difference  in  p(s) ,  p{s)  respectively,  at  different  values  of  their  arguments.  The 
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results  we  establish  in  the  sequel  make  use  of  energy  inequality  techniques  developed 
in  our  previous  work  [7,  8,  141.  In  particular,  in  [14]  (see  [14],  equation  (2.4)),  an 
inequality  of  the  form 

P\p{Q 2)  -  P(P2) I  <  c^piq2)^  -  p\ ,  (1.23) 

where  0  <  cx(p)  <  1 ,  was  used.  Using  a  similar  inequality,  we  are  led  to  assume 
here  that  the  first  pair  of  terms  on  the  right  in  (1.22),  multiplied  by  p ,  are  bounded 
by  the  first  term  on  the  right  in  (1.20).  The  second  two  terms  on  the  right  in  (1.22), 
multiplied  by  p ,  are  then  assumed  to  be  bounded  by  the  second  term  on  the  right 
in  (1.20),  with  e  (0  <  e  «  1)  being  a  measure  of  the  perturbation  in  p  generated 
by  p .  As  an  example,  suppose  that 

p(q2)  =  p(l+e2q2)-'12 ,  P{P2)  =  P ,  (1.24) 


so  that  u  satisfies  an  equation  of  minimal  surface  type  while  v  is  a  harmonic  func¬ 
tion.  The  constitutive  function  p  in  (1.24)  is  a  special  case  of  the  power-law  material 
(1.2)  with  n  =  1/2  and  b  =  e2/2.  Thus,  (1.24)  may  be  viewed  as  a  perturbation 
away  from  b  =  0  in  (1.2),  with  the  value  of  n  fixed  at  n  -  1/2 .  It  is  shown  in 
Appendix  A  that  (1.20)  is  satisfied  by  the  p ,  p  given  in  (1.24)  with 

Vl(p,e)  =  ep,  y2(p,e)=p2/ 2.  (1.25) 


In  Appendix  A  and  in  Sec.  3  it  is  also  shown  that  the  choices  for  yfp ,  s) ,  y2(p ,  e) 
given  by  (1.25)  do  indeed  satisfy  (1.18),  (1.19).  Another  example  we  shall  consider 
is  the  power-law  material  (1.2)  with  hardening  exponent  rt  =  1  +  e2 ,  so  that  the 
material  is  close  to  being  neo-Hookean.  Thus,  we  have 

P(42)  =  P  ( 1  +  -~-J )  *  p(p2)  =  p,  (1.26) 

V  !+«/ 

and  so  we  wish  to  compare  u  with  the  harmonic  function  v  satisfying  (1.1 3)-(  1.1 5). 
It  is  shown  in  Appendix  A  that  (1.20)  is  satisfied  in  this  case  with  the  choice 


7t(p,e)  =  — x{[e2/(l  +  s2)]'2/2  +Pb'/2} , 

1  -t-  e 


y2(p>£)  = 


bp2e2 
1  +  e2 


(1.27) 


In  Sec.  3  it  is  shown  that  (1.18),  (1.19)  again  hold  for  suitably  chosen  c,(e)  and  c2 . 

The  plan  of  the  paper  is  as  follows:  In  the  next  section,  we  establish  our  main 
result;  namely,  we  derive  an  exponential  decay  estimate  for  a  quadratic  functional 
defined  on  the  difference  between  u  and  v  .  We  show  that  this  “weighted  energy” 
is  bounded  above  by  an  exponential  that  decays  with  the  axial  variable  x,  and  we 
obtain  an  estimate  (lower  bound)  for  the  decay  rate.  Furthermore,  we  show  that 
the  energy  is  of  order  e4  ,  for  all  x, .  We  give  three  illustrative  examples  in  Sec. 
3,  two  of  which  involve  comparison  of  u  with  harmonic  functions,  while  the  third 
example  compares  solutions  for  two  softening  power-law  materials.  The  extension  to 
Neumann  boundary  conditions  is  described  in  Sec.  4. 
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2.  An  energy  comparison.  In  this  section,  we  establish  our  main  comparison  result 
between  u  and  v .  Let  w(x] ,  x2)  be  defined  by 

w{xx ,  x2)  =  u(xx ,  x2)  -  v{xx ,  x2).  (2.1) 

We  show  that  the  energy  measure 

E(z)  =  P{<i2)w_aw  a  dA,  (q2  =  u  pu  p) ,  (2.2) 

contained  in  the  subdomain 

R  2  =  {(JCi  ,  *2)|0  <2<x,<00,  0<JC2<fi} 

has  exponential  decay  in  z  and  is  of  order  e4  as  e  ->  0 .  In  fact,  we  shall  show  that 

E(z)<e4Cae~2k°z ,  z>  0,  (2.3) 

where  a  =  1  or  2  in  Case  1  or  2  respectively,  and  the  estimated  decay  rate  2ka  and 
the  amplitude  Ca  are  constants  which  can  be  explicitly  determined.  Note  that  the 
summation  convention  is  not  used  in  (2.3). 

The  result  (2.3)  is  established  in  several  stages.  First,  for  a  =  1  or  2  we  derive 
the  differential  inequality 

F'(z)  +  2dKa(z)F(z)<De4Ka{z)E(z),  z>  0,  (2.4) 

where  the  prime  denotes  differentiation  with  respect  to  z .  Here 

F(z)  =  E(z)  +  e4c2(l  +  cx)~'  E(z)/2,  (2.5) 

where  c,(e)  and  c2  are  the  quantities  introduced  in  (1.18),  (1.19)  respectively,  and 

£(*)=  /  P(P2)Vtav  adA,  (p2  =  vJ)vJj)  (2.6) 

jRz 

is  an  energy  measure  defined  on  solutions  v  of  the  comparison  problem  (1.12)- 
(1.15).  The  quantity  k  (z)  in  (2.4)  is 


Ba{z) 


(o=l,2). 


where 


1  Jo  ^ 


{q  )dx2  (Case  1), 

~'(q2)dx2  (Case  2). 


The  constants  d  and  D  in  (2.4)  are  given  by 


d=- Izh  , 

4(1  +c,)’ 


(3  +  c2)c22 


4(1 +C,)  4(1  +  c,)  (1  —  c,) 

respectively.  The  differential  inequality  (2.4)  is  integrated  once.  The  next  step  uses 
an  exponential  decay  estimate  for  E(z) ,  established  in  [7]  (see  (2.1)  of  [7]).  The 
estimate  is 

£(:)<c/‘!,  z  >0,  (2.10) 
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where 

and 


°a  =  ™anlMah  (a  =  1 ,  2) , 


Ga  =  £0exp 


2™cKaE o 


A/:A 


2 1.2 


(2.11) 


(2.12) 


In  (2.12)  E0  =  £(0)  denotes  the  total  energy  (assumed  finite)  contained  in  the  semi¬ 
infinite  strip  R .  Bounds  on  £(0)  in  terms  of  the  boundary  data  ( 1 . 1 5)  are  obtained 
in  [7].  When  (2. 1 0)  is  inserted  in  the  integrated  form  of  (2.4),  the  resulting  differential 
inequality  may  again  be  integrated.  This  introduces  the  total  weighted  energy  E0  = 
E(0)  which,  in  turn,  can  be  bounded  in  terms  of  E(0)  (see  (2.50)  below).  The  final 
step  makes  use  of  the  hypotheses  (1.7),  (1.8)  in  Cases  1,  2  respectively,  to  obtain  the 
desired  result  (2.3). 

To  establish  (2.4)  we  proceed  as  follows.  If  Lz  denotes  the  line  segment  jc,  =  z , 
0  <  x2  <  h ,  we  find,  by  using  the  divergence  theorem  and  (1.1),  (1.4),  (1.5),  (1.1 2)— 
(1.14)  that 


E(z)  —  ~  f  p(Q2)wwldx2  -  f  [p(q2)~  piP2)]wv  xdx2 
J  Lx  J  Lt 

-  JK  [p(q2)  -  P(p2)]v,fiW'fi  dA , 
where  we  use  the  notation 

y/dx2  =  ip(z,x2)dx2. 

Thus,  from  (2.13)  and  the  definition  of  p  in  (1.12),  we  have 

E{z)<  f  p{q2)\w\\w  x\dx2  +  f  \p(q2)- p(p2)\w\\v  t\dx2 
Jl ,  Jlz 

+  Ir  ~  P(P2)\P(W ,pw ,/?)1/2  dA  • 


(2.13) 


(2.14) 


(2.15) 


Denote  the  third  integral  on  the  right  in  (2. 1 5)  by  /3 .  Then  we  can  use  the  hypothesis 
(1.20),  together  with  (1.18)  and  (1.19)  to  obtain 


h  <  c,  Jk  \q  -  p\p(q2){w  fiw  j,)'12  dA 

+  c2e2  J  [p(q2)U2(w  pw  p)'/2p(p2)',2p]dA . 


By  virtue  of  the  definitions  of  p ,  q  ,  and  w  we  have 


2  ,  ,1/2,  ,1/2  ^ 
P  =  V,0U,0.  PI  =  (» (« 


(2.16) 


(2.17) 


and  so 


\q-P\  =  (q2-2pq+p2)'12  <  (v  pv  fi-2v  Bu  B+u  Bu  B)l/2  =  (w  Bw  B)l/2 .  (2.18) 


,r.r  .r.^ 


j 
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By  inserting  (2.18)  in  the  first  integral  on  the  right  in  (2.16),  using  Schwarz’s  inequal¬ 
ity  in  the  second  integral,  and  recalling  the  definitions  of  E(z) ,  E(z)  from  (2.2) 
and  (2.6)  respectively,  we  obtain 


/3  <  c}E(z)  +  c2e2El/2(z)El,2(z). 


(2.19) 


The  second  integral  on  the  right  in  (2.15),  denoted  by  l2 ,  is  bounded  in  a  similar 
fashion.  Thus,  since  |v  ,|  <  ( v  =  p ,  on  using  (1.20),  (1.18),  (1.19)  we  find 


72<^/  \q-p\p(q2)\w\dx2  +  c2e2  [  p(q2)'/2\w\p(p2)l/2p dx2 

Jl,  Jl,  (2.20) 

=  +  J2. 

Making  use  of  (2.18)  to  bound  7, ,  we  obtain 

Ji-cifL  p(q2)\w\(tv  pw  p)l/2dx2 

/  yn  /  \>/2  (2.2i) 

“C‘  \Jl  p^w-fiw<0dx2  J  [jL  p{Ql)w2  dx2j  . 

The  last  step  in  (2.21)  follows  from  Schwarz’s  inequality.  Using  a  scheme  involving 
a  change  of  variable  introduced  in  [7]  (see  pp.  314-315  of  (7])  (and  also  used  in  [8], 
[14]),  we  can  show  that 


(2.21) 


f  p(q2)w2 dx2  <  ^-1—  j  [P(q2))  'w22dx2  (Case  1), 
Jlz  7T  Jl,  ' 

f  \p{q2)]~lw2 dx2  <  f  p(q2)w22dx2  (Case  2), 
Jl,  n*  Jl, 


(2.22) 


where  Ba(z)  (a  =  1 , 2)  is  defined  in  (2.8).  By  using  (2.22)  and  the  left-hand  sides 
of  (1.7),  (1.8),  we  find  from  (2.21)  that 


J  <  C\^a,(Z)  f  p(q2)W  W  dx 

man  Jl,  •"  2 


(2.23) 


which,  in  view  of  (2.2)  and  (2.7),  can  be  written  as 

J.Sc^'izM-E'iz)).  (2.24) 

A  bound  for  J2 ,  defined  in  (2.20),  follows  similarly.  Thus,  by  using  Schwarz’s 
inequality,  we  have 


J2  <  c2e  .  p{q  )w  dx2 


p{p2)p2dx2 


(2.25) 


Again  using  (2.22),  the  left-hand  sides  of  (1.7),  (1.8),  the  definition  (2.6)  of  E ,  and 
the  obvious  inequality 


p(q2)w22dx2<  p{q2)w  ,w  „dx 2, 


(2.26) 
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we  find  from  (2.25)  that 

h  <  c/K-\z){[-E\z)]m[-E\z)\m} .  (2.27) 

Thus,  on  combining  (2.24),  (2.27),  we  obtain  from  (2.20)  that 


I2<clKa\z)[-E’(z)]  +  c/K;\z)[-E\z))'l2[-E'(z)]112.  (2.28) 

It  remains  to  obtain  an  upper  bound  for  /, ,  the  first  integral  on  the  right  in  (2. 1 5). 
By  using  Schwarz’s  inequality,  we  have 


(/t  p{ql)w2x  dx^j 


(2.29) 


Again  employing  (2.22),  the  left-hand  sides  of  ( 1 .7),  ( 1.8),  and  the  arithmetic-geomet¬ 
ric  mean  inequality,  we  find  that 

=  *:'(*) l-£'(*)]/2.  (2.30) 

Combining  the  results  (2.19),  (2.28),  and  (2.30)  on  the  right-hand  side  of  (2.15),  we 
find  that 

(1  -  Cl)E(z )  <  c2e2El,2(z)E(z)i/2  +  k~\z)(  1  +  2c,)[-£'(z)]/2 
+  c2e2K-a\z)[-E\z)]xl2[-E\z)\'12 . 

The  final  step  in  establishing  (2.4)  from  (2.31 )  is  carried  out  by  using  the  weighted 
arithmetic-geometric  mean  inequality 

lab  <  da 2  +  \b2  (8  >  0) 

o 

in  the  first  and  third  terms  on  the  right  in  (2.31).  Thus,  we  obtain 

(1  -c,  -  8l/2)E(z)  +  1/2(1  +  2c, +  S2)k~'(z)E'(z) 

+  c22e*K^(z)E'(z)/{232)-c22e*E(z)/(2Sl)  <  0 
for  arbitrary  <5, ,  S2  >  0 .  Choosing  <5,  =  1  -  c,  and  rearranging  (2.33)  we  get 


(2.32) 

(2.33) 


E'(z)  +  (1  +  2c,  +<52r'(l  -cl)Ka{z)E(z)  +  e4cl(  1  +  2c,  +  d2)~'s~' E'(z)  ^ 

-e*c22(l-ciy'(l  +  2ci+82)~lKa(z)E(z)<0 

for  arbitrary  S2  >  0 .  For  algebraic  simplicity  in  what  follows,  we  take  S2  =  1  in 
(2.34).  (Other  choices  for  ,  S2  may  also  be  made;  see  the  discussion  in  Sec.  3.) 
On  defining  F(z)  as  in  (2.5),  it  is  then  readily  verified  that  (2.34)  has  precisely  the 
desired  form  (2.4). 

The  differential  inequality  (2.4)  can  be  written  as 


F(z)expl^2dJ^  Kn(rj)dti}  <  De4Ka(z)E{z)exp  l^2d  j  Ka(rj)dr]^  ,  (2.35) 
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and  so,  on  integrating,  we  find  that 

F(z)  <  F{ 0) exp  {-2d  f  Ka(tj) dri  1 

!  <  r  t  <2-36> 

+  £>c  Ka(o)E(o)exp^-2d  J  Ka{r\)dq^  da . 

Now,  by  using  Schwarz’s  inequality,  we  can  easily  show  (cf.  pp.  315-316  of  [7])  that 

f  xa(n)dri  >  (z  -  a)2 j J'  Ba(s)ds  .  (2.37) 

Making  use  of  the  hypotheses  ( 1 .7)  or  ( 1 .8),  in  Cases  1  or  2  respectively,  we  can  show 
(see,  e.g.,  p.  316  of  [7])  that  (2.37)  leads  to 


where 


E(z)  =  f  p{q2)q2  dA  ( q 2  =  u  au  J 

Jr. 


(2.38) 


(2.39) 


is  the  energy  associated  with  the  problem  (1.1),  (1.4)-(  1.6).  Since  E(a)  is  decreasing 
in  a  (as  is  shown  in  [7]),  and  discarding  the  last  positive  term  in  (2.38),  we  obtain 


k  (q)dq  >  2  **  -  , 

man  Ja  aK  ~  Mah  mW 


(2.40) 


where  EQ  =  E( 0)  is  the  total  energy  (assumed  finite)  in  R .  By  inserting  (2.40)  (when 
a  =  0)  into  the  first-term  on  the  right  in  (2.36),  we  obtain 


F(z)  <  F( 0)exp 


dm  K  n 

_ a  n 

Mlh 2 


exp 


-2  dmnnz 

Wh 


De*  Ka(o)E{o)exp{-2d  J  Kn(q)dq{  da. 


(2.41 ; 


The  exponentially  decaying  term  in  the  first  expression  on  the  right  in  (2.41)  has 
a  decay  rate  2va  (a  =  1,2)  given  by 


v  = 


dm  n 

Hth' 


(2.42) 


where  d  is  given  in  (2.9), .  We  now  use  the  decay  estimate  for  E(a) ,  given  by  (2.10), 
in  the  integrand  in  the  integral  in  (2.41).  This  substitution  will  then  introduce  the 
decay  rate  ua  ,  given  by  (2.1 1),  into  the  final  estimate.  By  using  integration  by  parts, 
we  find  that  the  entire  second  expression  on  the  right  in  (2.41),  denoted  by  J ,  is 
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such  that 


J  < 


DGe 4 


-  Id 


3^-{  e  ^"exp  j-2 d  J  *,(>?)  d>?J 


+  2i>alo  6  2V7eXp{~2rf  /  <**} 

exp  {-2rf  jf  *,(,  )*}] 


-  2d 
+ 


(2.43) 


DG  e4i> 


r  e  v  [z  -io  , 
_1> _ a  /  o  L° 

Id  J0 


exp 


dmaKnE,  \ 

mW 


exp 


-2  dm(z-a) 


M  h 


da , 


where  we  have  used  (2.40)  to  obtain  the  final  inequality  in  (2.43).  By  discarding  the 
second  (negative)  term  in  (2.43)  and  inserting  the  resulting  bound  in  (2.41),  we  find 
that 


F(z)<F(  0)exp 


\2vK„\ 

Q  Q 

M  h 

a  J 


-2v  :  t  -2t  z 

e  "  H — ~v" — e 


2d 

-2^2 


.  JgA 

2d  CXP  MqA  P0-*0 


(2.44) 


where  if  va  =  va,  then  the  limit  is  understood  in  the  last  term  of  (2.44). 

The  desired  decay  estimate  (2.3)  follows  from  (2.44)  by  observing  that  the  ex¬ 
ponentially  decaying  terms  on  the  right-hand  side  of  (2.44)  can  be  bounded  by  an 
exponentially  decaying  term  of  the  form  exp(-2 kaz) ,  where 

K  =  minK  ,*«)  =  £  min  •  (2.45) 


We  recall  from  (2.9),  that  d  =  (1  -c,)/[4(l  +c,)]  and  c,  =  c,(e)  is  the  quantity 
defined  by  (1.18),  (1.20).  To  obtain  the  sharpest  possible  estimated  decay  rate  ka  , 
one  should  choose  ma ,  ma  as  large  as  possible  and  Ma  ,  Ma  as  small  as  possible. 
In  the  illustrative  examples  discussed  in  Sec.  3,  we  have 


va<va.  (2.46) 

In  fact,  it  is  shown  in  Appendix  C  that  (2.46)  holds  if  the  constitutive  functions  p 
and  p  satisfy  a  condition  suggested  by  (1.20)  (see  equation  (C.l)).  On  assuming  that 
(2.46)  holds,  and  observing  from  (2.5)  that  F{z)  >  E(z) ,  we  find  from  (2.44)  that 


where 


E(z)  <  Hae 


-Ivz 


s4DG.. 


z  >  0, 

(2.47) 

u  ( 2v  K  Eu\ 

"Jexpv  M"h  )■ 

(2.48) 

Observe  that  the  quantity  Ga  defined  in  (2.12)  contains  E0 ,  and  that  F( 0) ,  given 
by  (2.5)  with  z  =  0 ,  contains  both  E0  and  E0 .  We  can  easily  obtain  an  upper  bound 
for  E0  in  terms  of  £0.  By  virtue  of  the  definition  (2.1)  of  w  and  by  using  (1.6), 
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(1.15),  we  see  that  to  =  0  on  x,  =  0.  Thus,  (2.15)  evaluated  at  z  -  0  yields 
E0  <  I3(0)  and  so,  by  (2.19),  we  find  that 


(1  ~c,)E0<  c2t 


£l/2£1/2 
c0  ’ 


(2.49) 


and  so 


£0< 


2  4 

c2e 


(1  -ct) 


E 

2  c0 ' 


(2.50) 


Thus,  we  have  from  (2.5)  that 

F( 0)  =  E0  +  c2V(l  +c,rI£0/2  <  c2Y(3  +  cf)(l  +  c,f '(1  -  cx)~2EJ 2.  (2.51) 

Upper  bounds  for  the  total  energies  E0  and  E0  in  terms  of  the  boundary  data  have 
been  established  in  [7],  and  these  bounds  can  be  used  in  (2.50),  (2.51),  and  (2.48)  to 
obtain,  from  (2.47),  the  desired  result 

E(z)  <  tCae~1VaZ ,  z  >  0  (2.52) 


for  a  computable  constant  Ca . 

If  the  hypothesis  (C.l)  of  Appendix  C  is  not  made  (so  that  (2.46)  is  no  longer 
assured)  one  can  still  proceed  from  (2.44),  and  use  (2.50),  (2.5 1 )  to  obtain  an  estimate 
of  the  form 

E(z)  <  e*Cae~2k°z ,  z  >  0  (2.53) 

for  a  computable  constant  Ca ,  where  ka  is  given  by  (2.45).  This  completes  the 
derivation  of  (2.3). 


3.  Illustrative  examples.  In  this  section  we  illustrate  by  three  examples  how  to 
make  explicit  the  comparison  (2.52)  of  the  solution  of  problem  ( 1 . 1 ),  ( 1 ,4)-(  1 .6)  with 
the  solution  of  the  perturbed  problem  (1.12),  (1.1 3)— ( 1 . 15).  Two  of  these  examples 
were  mentioned  in  Sec.  1,  i.e.,  the  p  and  p  given  by  (1.24)  and  (1.26).  We  consider 
first  the  problem  defined  by  (1.24). 

Problem  1.  Here  p{q2)  and  p(p2)  are  given  by  (1.24),  i.e., 

p(q2)  =  p{l  +  e2q2)~l/2 ,  p{p2)  =  p.  (3.1) 

We  observe,  in  connection  with  (3.1),  that  one  might  think  of  p  as  a  perturbation  of 
p  rather  than  p  being  a  perturbation  of  p ,  but  it  is  immaterial  which  we  designate 
as  the  perturbation  of  the  other.  At  any  rate,  as  mentioned  in  Sec.  1,  it  is  shown  in 
Appendix  A  that  (1.20)  is  satisfied  with  y{{p,  e)  and  y2(p ,  e)  given  by  (1.25),  i.e., 

yt(p,e)  =  ep,  y2(p ,  e)  =  p2 / 2.  (3.2) 

In  order  to  find  the  constants  c,  (e)  and  c2  defined  by  ( 1 . 1 8)  and  (1.19)  respectively, 
it  is  necessary  to  find  a  bound  for  the  maximum  value  of  p  =  |gradi>|  on  7t .  Since 
p  is  the  gradient  of  a  harmonic  function  this  is  not  difficult  to  accomplish  with  the 
conditions  on  f{x2)  implied  by  the  fact  that  v  e  C2{R )  n  C1  ( R ) .  The  derivation  of 
an  appropriate  bound  is  carried  out  in  Appendix  B  where  it  is  shown  that,  with  the 
notation 

max p2  =  P2 , 

R 


(3.3) 
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then 

+  (34) 

It  follows  that  cx  (e)  and  c2  can  be  taken  as 

cl(e)  =  eP,  c2  =  P2/ 2.  (3.5) 

Thus,  for  sufficiently  small  e  ,  the  last  inequality  in  (1.18)  is  satisfied. 

The  p  and  p  given  by  (3.1)  belong  to  Case  2  (see  (1.8)  and  (1.17)).  Clearly  (1.17) 
is  satisfied  with  the  choices 


m2  =  p 


-l 


M2  =  p 


Since 


/.  .  2  2,1/2  1  +  e2q2 

(l+eq)  =  ■■■■  ?  ->.,/•>  $  1  + 


o 

II 

(3.6) 

2  2 

E  q 

(1  +eV],/2’ 

(3.7) 

2  -2 
,  =  e  P  , 

(3-8) 

it  follows  that  (1.8)  is  satisfied  for  the  choices 

. — i  w  -i 


(cf.  [7],  p.  325).  Recalling  the  definitions  of  02 ,  u2  from  (2.1 1),  (2.42)  respectively, 
we  have 

*  *d  (3.9) 


V2  ~  h  ’ 


h 


with  d  given,  by  virtue  of  (2.9)  and  (3.5),  by 

(1-eP) 


d  = 


(3.10) 


4(1  +eP) 

Thus,  it  is  clear  that  v2  <  t>2  so  that  (2.46)  holds  and  the  decay  estimate  (2.47),  with 
a  =  2 ,  holds  for  Problem  1  with 


i/,  = 


n  (1  -eP) 


2  4h(\+eP)' 

Since  K2  =  0 ,  the  expression  (2.12),  with  a  =  2 ,  simplifies  to 


G2-  E0. 


(3.11) 


(3.12) 


Also,  since  by  (3.8)  the  constant  K2  is  of  order  e2 ,  the  decay  estima*e  (2.47),  (2.48) 
can  be  written  as 

E{z)  <e4LE0[l+O{e2)]e~2,/2Z  (3.13) 

where  v2  is  given  by  (3.1 1).  Here  the  constant  L  is  given  by 


2  L  = 


n*» 

^-(3  +  P2e2)(l  +  Pe)-I(l  -  Pef 2  +  D/(d(  1  -  d)) 


(3.14) 


where  D  is  defined  by  (2.9).  For  Problem  1,  on  using  (3.5),  (3.6),  (3.8),  and  (3.10), 
we  have 

(1-eP) 


D  = 


(3  +  e2P2)P4 


16(1  +eP)z(l  ~eP) 


d  = 


4(1  +eP)  ’ 


(3.15) 
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We  note  that  the  decay  rate  given  by  (3.1 1)  is  much  too  conservative  since  a  simple 
estimate  based  on  the  arithmetic-geometric  mean  inequality  gives 


E(z) 


<2  |£(z)  +  piq^v  av  adA^ 
<2{E(z)  +  E(z)} 

-2xz/h 


(3.16) 


<  Ne 


The  last  inequality  follows  from  results  established  in  [7],  Here  N  is  a  constant 
which  is  0(1)  in  the  parameter  e .  We  should  point  out,  however,  that  in  this  paper 
we  are  emphasizing  the  feature  that  E(z)  is  0{eA)  for  all  z  and  hence  that  u  and 
v  are  close  in  energy  measure  for  all  z  .  At  various  steps  in  our  proof  we  have  used 
inequalities  that  preserve  the  0{e4) ,  but  which  result  in  a  poorer  estimate  of  the 
actual  decay  rate.  For  instance,  in  this  problem  c,(e)  is  O(e) .  If  we  were  interested 
in  improving  the  estimate  of  the  decay  rate,  then  in  (2.33)  we  could  choose 


—  *  ^2  —  &2E 

for  arbitrary  positive  constants  /?,  and  /?2  .  Then  (2.33)  would  become 
[\-(P  +  fii/2)e)E{z)  +  (\/2)[\+{2P  +  P2)e)K-\z)E'{z) 

+  P*M2)-'z\-\z)E\z)  -  PA(&fiirleiE(z)  <  0. 

By  continuing  the  proof  as  above  we  would  obtain,  instead  of  (3.13), 

E(z)  <£3C2V2<Cj*2, 

where  C2  is  a  computable  constant,  but  now 

.  *  f  \-e[P  +  0l/2]\ 

2  h\  l+e[2P  +  02)  /• 

Thus  by  sacrificing  a  power  of  e  in  the  multiplicative  factor  in  the  decay  estimate  we 
are  able  to  improve  the  estimated  decay  rate  obtaining  a  new  estimate  which  differs 
from  the  expected  decay  rate  by  a  term  of  0(e) .  This  trade  off  between  the  size 
of  the  multiplicative  factor  in  the  estimate  and  the  decay  rate  is  reminiscent  of  that 
observed  by  Horgan  and  Payne  [7,  8]  in  related  contexts. 


Before  proceeding  to  the  next  example,  we  give  an  indication  of  how  Problem  1 
might  arise  in  applications.  Let  <t>  be  the  solution  of  the  following  minimal  surface 
problem: 


{U  +  |V<£|2r'/20a}.a  =  O  on  R 

(3.21) 

with  the  boundary  conditions 

</>(*,,  0)  =  0,  <f>{xl,h)  =  0,  x,  >  0, 

(3.22) 

0,0  |  — ►  0  (uniformly  in  x2)  as  jc,  -►  oo , 

(3.23) 

0(0,  x2)  =  ef(x2) ,  0  <  x2  <  h . 

(3.24) 

Then  if  we  set 

0  =  e«. 

(3.25) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 
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the  function  u  will  be  the  solution  of  (1.1),  (1.4)— (1.6)  with  the  p(q2)  of  (1.24), 
the  case  we  have  just  considered.  Thus,  we  might  think  of  approximating  by  ev 
where  v  is  the  harmonic  function  taking  the  value  zero  on  the  long  sides  and  the 
value  f(x2)  on  the  end  x,  =  0.  A  measure  of  the  merit  of  this  approximation 
follows  from  the  observation  that 


E\z)=  /  [l  +  \V4>)2)-l,2\V(4>-ev)\2dA 

J*, 

=  e2E(z) 

<  eC2e~2v'z ,  z  >  0. 

In  the  last  step  we  make  use  of  (2.52)  with  a  -  2. 

Problem  2.  Here  p  and  p  are  given  by  (1.26),  i.e.. 


(3.26) 


p{q2)  =  p[l  +  y£2\  .  P{P2)  =  P-  (3.27) 

As  we  mentioned  in  Sec.  1 ,  the  problem  for  u  may  be  thought  of  as  a  problem  for 
a  material  that  is  close  to  neo-Hookean  while  the  problem  for  v  is  the  neighboring 
neo-Hookean  problem  governed  by  Laplace’s  equation.  As  we  observed  in  Sec.  1,  y, 
and  y2  are  now  given  by  (see  Appendix  A) 


h(P,  «)  =  ^^{[e2/(l  +  e2)f'2+pb'12}, 
1  +  e 


V2(P>£)  = 


bp2e2 

l+e2 


(3.28) 


Since  the  v  in  the  perturbed  problem  is  again  a  harmonic  function  we  can  use  (3.4) 
as  the  bound  for  P ,  the  maximum  value  of  p  in  R .  We  thus  arrive  at 


.  ,  b  e  .  „  ,i/2 

(fi)  —  2^"^^  E  }  , 

l+e 


C2  = 


bP2e 2 
1  +  e2 


(3.29) 


Thus,  for  sufficiently  small  e ,  the  last  inequality  in  (1.18)  is  satisfied. 

The  p  and  p  given  by  (3.27)  belong  to  Case  1  (see  (1.7)  and  (1.16)).  Clearly 
(1.16)  is  satisfied  with  the  choices 


mx=p,  Mx=p,  Kx=0. 


Since 


1  < 


1  + 


bq 2 


d+^) 


<  1  + 


it  follows  that  (1.7)  is  satisfied  for  the  choices 


bq2e2 

l+e2 


be2 


mx=p,  Mx  =p,  Kx=  s 

1  +  el 


(cf.  [7],  p.  324).  Thus,  we  have 


v\  =  T. 


n 

h 


nd 


1,1  =  h  ’ 


(3.30) 

(3.31) 


(3.32) 


(3.33) 
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where  d  is  given  by  (2.9)  and  (3.29), .  Again  we  see  that  vx  <  so  that  (2.46) 
holds  and  so  the  decay  estimate  (2.52)  (with  a  =  1)  holds  for  Problem  2  with  ux 
given  by  (3.33)2 .  Since  A?,  =  0  and  K{  =  0(e2) ,  we  can  again  obtain  an  estimate  of 
the  form  (3.13),  with  v2  replaced  by  vx  for  Problem  2.  Also,  since  c, (e)  is  0(e2)  , 
it  is  again  possible  to  improve  the  estimated  decay  rate  at  the  expense  of  lowering 
the  order  of  e  in  the  multiplicative  factor  by  making  the  appropriate  choices  for  <5, 
and  S2  in  (2.33). 

Problem  3.  As  a  third  example  we  consider  the  following  expressions  for  p  and 


P- 


p(q2)  =  p 


p(p2)  =  p[i  +  p2)  1/4  • 


(3.34) 


We  observe  that  these  values  of  p  and  p  are  of  the  power-law  form  (1.2)  with 


b  = 


3 

4(1  +e2) ' 


(3.35) 


Thus,  we  are  considering  softening  materials  with  a  fixed  value  of  the  softening 
parameter,  and  the  perturbation  may  be  viewed  as  a  perturbation  in  the  material 
parameter  b  away  from  the  value  3/4.  We  observe  that  this  problem  differs  from 
the  previous  examples  in  that  the  differential  equation  for  v  is  now  also  quasilinear. 

It  is  shown  in  Appendix  A  that  for  this  example 


y,(p,e)  = 


(iWr,/4/> 

(1  +p2)l/2  +  [(1  +  e2)(I  +p2)]l/4 


y2(p,e)  =  (l+£2rV/4.  (3.36) 


To  find  the  constants  c,(e)  and  c2  ,  defined  by  (1.18),  (1.19)  respectively,  one  would 
require  a  bound  for  P,  the  maximum  of  p  =  |  grad  v  |  on  R .  In  contrast  with  the 
analogous  situation  for  Problems  1  and  2  where  such  a  bound  was  obtained  in  (3.3), 
(3.4),  here  one  requires  an  estimate  for  the  gradient  of  solutions  of  the  quasilinear 
differential  equation  (1.12),  with  p  given  by  (3.34)2 .  While  techniques  for  obtaining 
such  bounds  are  known,  the  detailed  calculations  required  to  obtain  explicit  estimates 
are  beyond  the  scope  of  the  present  work. 

From  (3.36),  it  follows  that  c,(e)  and  c2  may  be  taken  as 


c,(e) 


_ (1  +e2)~1/4P _ 

(1  +  P2),/2  +  [(1  +  e2)(l  +  />2)),/4 


c2  =  (l+e2)  V/4.  (3.37) 


For  e  sufficiently  small,  the  last  inequality  in  (1.18)  is  satisfied. 

Since 

,  2. —  1/4, ,  ,  2  ,  2.1/4 

p  ( q)  =  p  (1+e)  '(*+«+«) 

=  /T'(l  +e2)"'/4(l  +q2+e2)(l  +q2  +  e2)~3'4  (3.38) 

=  /*_l(l  +  e2)3/4(l  +  Q2  +  £2r3/4+/T2(  1  +  e2)-l/2(l  +  ,2  +  e2)~l,2q2p(q2). 


the  conditions  (1.8)  of  Case  2  hold  with 

m2  =  p-',  M2  =  p~',  K2  =  p~2(\ +e2)~' 


(3.39) 
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Also,  as  was  shown  in  [7]  (see  pp.  324-325  of  [7])  the  conditions  (1.17)  of  Case  2 
hold  with 


(3.40) 


Thus,  we  have 


with  d  given  by 


d  = 


1 

4(1  +c,r 


(3.41) 

(3.42) 


and  c,  given  by  (3.37), .  Again,  we  have  v2  <  02  so  that  (2.46)  holds  and  the  decay 
estimate  (2.52),  with  a  =  2,  holds  for  Problem  3  with  v2  given  by  (3.41),  (3.42). 

The  decay  behavior  in  this  problem  is  somewhat  different  from  that  in  the  two 
previous  examples.  In  the  first  two  problems  p  -*  p  as  q2  -*  0  whereas  in  this 
example  that  feature  is  no  longer  true.  We  observe  also  that  since  c,(c)  is  0(1) 
in  e  as  e  -*  0  we  no  longer  have  the  possibility  of  a  trade  off  between  the  decay 
coefficient  and  the  order  of  e  in  the  multiplicative  factor.  At  the  same  time,  since 
k2  is  different  from  zero,  the  constant  C2  in  (2.52)  will  in  general  be  much  larger 
than  in  the  first  two  examples. 

Finally,  we  provide  another  interpretation  for  the  perturbation  involved  in  Prob¬ 
lem  3.  Let  0  be  the  solution  of  the  problem 


{[1  +  |V0|2]  ,/4fa},o  =  0  on  R  (3.43) 


subject  to  the  boundary  conditions 


0(xl,O)  =  O,  <t>(xl,h)  =  0,  *,>0,  (3.44) 

0 ,  0  ,  0  (uniformly  in  jt2)  as  Jt,  -*  oo ,  (3.45) 

0(0,  jt2)  =  (1  +  e2)~l,2f(x2),  0  <  x2<  h .  (3.46) 

Then,  if  we  set 

0  =  (l+e2fl/2n,  (3.47) 

u  will  be  the  solution  of  (1.1),  ( 1 .4)— ( 1 .6)  with  the  p  of  (3.34), .  Thus  we  might 
think  of  approximating  0  by  (l+e2)-l/2t;  where  v  satisfies  (1.1 1)-(  1.1 4)  with  p(p2) 
given  by  (3.34)2 .  A  measure  of  the  merit  of  this  approximation  can  be  assessed  from 
the  fact  that 


E**(z)=  f  [1  +  |V0|2]  ',4\V[<t>-(\+e2)-'l2v]\2dA  =  (l+e2)-'E(; 
Jr, 


<  e  r  p~2v iZ 

£  TT?  2  ■ 


(3.48) 


z  >  0. 


In  the  last  step  we  use  (2.52)  with  a  =  2.  Note  that,  in  contrast  to  the  result  (3.26), 
the  right-hand  side  of  (3.48)  is  of  order  e4  as  e  -»  0  rather  than  e6 ,  as  in  (3.26). 


458 


C.  O.  MORGAN  and  L.  E.  PAYNE 


4.  Concluding  remarks.  The  preceding  results  are  also  valid  when  the  Dirichlet 
boundary  conditions  (1.4),  (1.6)  are  replaced  by  appropriate  Neumann  boundary 
conditions  that  correspond  to  traction  boundary  conditions  in  the  antiplane  shear 
context.  Thus  we  now  consider  solutions  of  (1.1)  on  R  subject  to  the  boundary 
conditions 


«i2(x, ,  0)  =  0,  w2(x,,A)  =  0,  x,>0,  (4.1) 

u ,  u  |  *  0  (uniformly  in  x2)  as  x,  -»  oc ,  (4.2) 

p(q2)u  ,  =  g(x2)  onx,=0,  0  <  x2  <  A .  (4.3) 

The  prescribed  function  g ,  assumed  sufficiently  smooth,  must  also  satisfy  the  “self¬ 
equilibration”  condition 


[  g{x2)dx2  =  0, 

Jo 


(4.4) 


which  follows  from  application  of  the  divergence  theorem  to  (1.1)  and  by  using 
(4.1)— (4.3).  In  fact  the  same  argument  shows  that 


f  p(q2)u  t  dx2  =  0,  z>  0,  (4.5) 

on  every  cross  section  Lz  .  The  problem  for  v  is  again  governed  by  ( 1 . 1 2),  but  now 
subject  to  the  boundary  conditions 


v  2(x, ,  0)  =  0,  V  2(x, ,  A)  =  0,  X|>0, 

(4.6) 

v  ,v  ,  -♦  0  (uniformly  in  x2)  as  x,  —  oo , 

(4.7) 

$(p2)v  ,  =  g(x2)  on  x,  =  0,  0  <  x2  <  A . 

(4.8) 

As  in  the  problem  for  u ,  we  find  that 

j  p(p2)v  ,  dx2  =  0  z>  0. 

(4.9) 

By  combining  (4.5)  and  (4.9),  we  see  that 

f  [/>(<?2)«,1  ~  p{p2)v ,\]dx2  =  0,  z  >0. 

(4.10) 

L, 


The  proof  of  (2.3)  given  in  Sec.  2  remains  valid  for  the  Neumann  boundary  conditions 


up  to  equation  (2.13).  We  now  introduce  the  functions  IT^x, ,  x2)  (n  = 
Cases  1,  2  respectively)  by 

Wa  =  w-wa(xx). 

1 ,  2,  in 

(4.11) 

where 

1  fH  2 

«M*i)  =  J0  pi<q  ^wdxi' 

w2=  IT  f  [P(Q2)V'wdx2- 
"2  JO 

(4.12) 

By  virtue  of  their  definitions 

/  P(Q2Widx2  =  0, 

Jo 

[\p(q2)]''w2dx2  =  0. 

Jo 

(4.13) 
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By  using  (4.10),  we  can  see  that  (2.13)  can  be  written  as 

E(z)=  -  J  p(q2)Waw  tdx2-  J  [p(q2)~  p(p2)]Wav  tdx2 

-  j  I P{q2)  -  P(p2)]v  ew  fidA . 


(4-14) 


Thus,  for  the  Neumann  boundary  conditions,  (4.14)  can  be  obtained  from  (2.13) 
by  formally  replacing  w  in  the  line  integral  terms  by  Wa  (a  =  1 , 2  in  Cases  1 ,  2 
respectively).  One  now  proceeds  from  (4.14)  exactly  as  was  done  in  Sec.  2,  except 
that  in  (2.22),  the  w 2  terms  on  the  left  are  replaced  by  W2  (a  =  1 , 2) .  In  deriving 
this  version  of  (2.22),  one  makes  use  of  the  zero  average  conditions  (4. 1 3).  (A  similar 
procedure  was  followed  in  (7,  8].)  In  this  way,  one  arrives  at  the  estimates  (2.47), 
(2.48),  where  the  total  energies  E0,  E0 ,  and  E0  are  now  understood  to  be  those 
associated  with  the  Neumann  boundary  conditions.  With  this  understanding,  one 
finds  that  (2.49)-(2.53)  are  also  valid  for  Neumann  boundary-value  problems.  The 
techniques  described  in  [7,  8]  may  be  used  to  obtain  the  required  bounds  for  the  total 
energies  E0 ,  E0  in  terms  of  the  boundary  data  (4.3),  (4.8). 

Appendix  A.  Verification  of  (1.25),  (1.27),  and  (3.36).  If  p(q2)  and  p(p2)  are 
defined  as  in  (1.24)  we  observe  that 

p\p{q2)  -  p(p2)\  =  PPlH  +eVr1/2-  1] 

2  2 

_ _ ppe  q _ 

(1  +eV)l/2[l  +(1  +fiV)l/2]  (A-*) 

_  pe2[<?(g  -P)+  P(q  -p)+  P2 1 

A\l+eV),/2ll+(i+«V)l/2l' 


Clearly  then 


„  pep(q*)[eq  +  ep]  .  e2p2lp(q2)p(p2))'12  ^ 

^  )-pU’  "  £  (!  +(l  +«V)'/21  (I  +«J9!),,4[l  +d  +«V)1,J1 '  <A'2) 


Now  setting  eq  =  s  and  tp  ~t  we  seek  an  upper  bound  for 

t  \  (•$  +  0  /at 

ip(s)  =  1 - rjn  ■  (A.3 

1  +  (1  +s2)'/2 

For  /  <  1  it  is  easily  checked  that  tp  is  an  increasing  function  of  s  and  thus  that 


y/(5)  <  1  . 


This  leads  to 


p\p{q2)  -  p{p2)\  <  pep{q2)\q  -p\  +  —^—[p(q2)P(p2)]l/2  ■ 

By  comparing  with  (1.20),  we  observe  that 

yl(p,e)=pe,  y2(p ,  e)  =  p2 /2, 
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which  is  (1.25).  We  observe  further  that  ep  <  1  is  required  by  condition  (1.18). 
For  p  bounded,  this  merely  imposes  a  restriction  on  the  allowable  magnitude  of  the 
parameter  e  (see  the  discussion  in  Sec.  3). 

Next  we  establish  (1.27)  for  p(g2)  and  p(p2)  given  by  (1.26).  Now 

p\p(q2)  -  p{p2)\  =  /^|^i  +  ~l  >  (A-7) 

and  we  bound  the  expression  on  the  right  making  use  of  an  inequality  from 
{16,  p.  356]  (stated  in  a  somewhat  different  form),  i.e.,  for  any  real  number  a  and 
0<a<  1/2, 

a-  1  <a(a-  l)a~(,_")/2  (A.8) 

(cf.  (A.  12),  p.  131  of  [8]).  Thus 

p\p(q2)  -P(P2)\  <  /ibe2p2q2(l+e2r'[l+b2/(l+e2))-{'-'')12 

=  Pbe2p[(q  -  p)(q  +  p)  +  p2]  (A-9) 

(1  +  e2)[l  +  bq2/{l  4-  e2)](l_*2)/2 

It  follows  then  that 


(1  +  r) 


p  +  q _ 

+  bq2/(  l+e2)](,+£2)/2. 


,  bpV  [p(q2)p{p2)]'12 
1  +  e2  [1  +  b2q2/(l  +  e2)] 


(A.  10) 


(A.  15) 
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Thus,  we  find  from  (A.  12)  that  y,(p,  e)  can  be  chosen  as 

yl(p,e)  =  pb'/2e2/(l  +  e2){[e2/(l  +e2)]* 12  +pb'12}. 

Equations  (A.  16)  and  (A.  1 1  )2  give  precisely  the  expressions  (1.27). 
To  establish  (3.36)  we  start  with 

l(l+?J  +  e2),/4-(I+p2),/4(I+e2),/4| 


(A.16) 


p!p(<?2)  -  £(p2)  I  =  pp! 


(i  +q2  +  «2)l/4(i  +p2),/4 


=  p 


p|(l  +q2  +  e2),/2  -  (I  +  p2)l/2(l  +  e2)'/2| 


(1  +?2  +  *2)l/4(l  +p2),/4[(l  +92  +  e2),/4  +  (l+p2),/4(l  +«2)l/4] 


(A.  17) 


Further  simplication  leads  to 

pIp(«J)  -  /HP1) I 

_  _ PP\q2  -p2  -eVl _ _ 

(1  +  «2  +  «2),/4(l  +p2),/4l(l  +  «2  +  *2),/4  +  (l  +  P2),/4(1  +«J),/4)[(1  +  </2  +  *2)‘/2  +  (I  +  p2)1/20  +e2),/2l 

-(l+«V‘/W)l«-Pl 

x  _ P(«  +  P) _ 

(1  +p2)l/4[(l  +«2  +  e2)'/4  +  (l  +p2)1/4(I  +«2),/4l[(l  +«2+t2),/2  +  (l  +  p2),/2(l  +«2)‘/2] 

+ _ (H-e2)-,/2pV[p(g2)PO>2)}'/2 _ 

[(1  +  *2  +  «2)(t  +  P2)1,/,I<1  +  ?2  +  e2)l/4  +  (I  +  p2),/4(l  +  «2)1/4][( I  +  q2  +  «2)1/2  +  (1  +  p2)l/2( I  + 12)'/2]  ’ 

(A.  18) 


By  comparing  with  (1.20),  we  see  that  we  may  choose 

y,(p, «) 


=  max 


(l+e2r'/4p(<7  +  p) 


«  (1  +P2),/4[(1  +?2  +  e2)1/4+(l  +P2),/4(1  +  e2),/4][(l  +?2  +  e2),/2+(t  +p2),/2(I  +e2),/2J 


<  (1  +e2)_,/4 


max 


(1+P2)1/4  l  «  [(l  +  92  +  e2),/4+(l+P2),/4lJ  l  «  (l+(?2)l/2  +  (l+p2) 


}{ 


max 


p  +  q 


X*/2 


(l+e2)~l/4p 

"  (l+p2),/2  +  [(I+e2)(l+p2)J,/4’ 

where  the  second  quantity  in  braces  has  been  bounded  above  by  unity. 
Clearly,  also 

y2{p,e)  =  (1  +e2)_1/V 

x  max{[(  1  +  q  +  e2)(  1  +  p 2)]~'/S 

<1 

x[(l+,2  +  e2r  +  (l+p2)^(l+e2)''r 

xi(i+92+g2),/2+(i+p2),/2(i+£2)l/2r1} 

<(i±iVV 

4(1  +e2)3/4 
=  (l+e2rV/4. 

Inequalities  (A.  19)  and  (A.20)  yield  the  desired  expressions  (3.36). 


(A.19) 


(A.20) 


J 
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Appendix  B.  Proof  of  (3.4).  If  v  is  a  harmonic  function  that  vanishes  on  the  lateral 
sides  of  R  and  vanishes  at  infinity,  then  it  follows  from  HopPs  second  principle  that 
p2  -  |Vt;|2  must  take  its  maximum  value  on  the  end  z  =  0.  Furthermore,  Payne 
and  Philippin  [17]  have  shown  that  if  g  is  any  other  harmonic  function  in  R ,  then 
the  quantity  eag )Vv|2 ,  for  any  real  constant  a ,  also  takes  its  maximum  value  either 
at  x,  =  oo  or  at  some  point  on  the  boundary  of  R .  We  apply  this  result  making  the 
special  choice  g  =  2x, ,  that  is,  we  define 

9  =  e2ax'p2,  0  <a<£.  (B.l) 

n 


We  remark  that  since  p2  is  0(e~2xx'/h)  as  x,  -*  oo,  for  our  choice  of  a  we  are 
assured  that  the  maximum  does  not  occur  at  infinity.  It  is  also  easily  seen  that  the 
maximum  cannot  occur  on  x2  =  0  or  x2  =  h  .  To  see  this  suppose  that  the  maximum 
did  occur  at  a  point  Q  on  x2  =  h  ,  x,  >  0 .  Then  since  6^0  in  R ,  it  follows  from 
Hopfs  second  principle  that 

(B'2) 

But,  at  each  point  on  x2  =  h , 


d6 


~~  =  2(v  xv  n  +  v  2v =  2(t>  (v  I2  -  v  2v  n)e"ax>  =  0, 


2ax,  _ , 


(B.3) 


since  both  v  ,  and  v  ,,  vanish  on  x2  =  h.  Thus,  9  cannot  take  its  maximum 
value  on  jc2  =  h ,  jc,  >  0 .  A  similar  argument  shows  that  the  maximum  value  is 
not  taken  on  x2  =  0,  x,  >0.  Thus  the  maximum  must  occur  on  x,  =  0.  If 
the  maximum  occurs  at  (0,  0)  or  (0,  h) ,  then  the  assumed  continuity  (recall  that 
v  €  C2(R)  n  Cl(R))  assures  us  that 


p2  <rmx{[f(0)]2,[f'm2}.  (B.4) 

On  the  other  hand,  suppose  that  the  maximum  value  of  6  occurs  at  a  point  Qx  = 
(0,  x2)  for  0  <  x2  <  h.  Then  at  Q,  ,  we  have  (again  by  using  HopPs  second 
principle) 


vtvu+v2vi  2  +  aP2<0>  (B.5) 

V,lV,l2  +  V,2V,22  = 

By  using  the  differential  equation  and  the  boundary  condition  (1.15)  at  x,  =  0,  we 
may  rewrite  these  expressions  as 


-v  ,/"  +f'v  l2  +  ap2  <  0,  (B.7) 

v,lv,n  +  f'f"  =  ° •  (B.8) 

If  v  X{QX)  =  0,  then  by  (B.8)  it  follows  that  f"{Qx )  =  0  (since  f'(Qx)  /  0, 
otherwise  p  would  vanish  at  this  point).  This  would  again  lead  to 

/>2<[/'(0,)]2.  (B.9) 


CONSTITUTIVE  LAW  PERTURBATIONS 


463 


However,  if  v  ,(<?,)  ^  0,  we  may  eliminate  v  l2  between  (B.7)  and  (B.8)  to  find 
that  at  Qx 

p2{a  -  f  " /v  x)  <  0  (B.  10) 

or 

v\(Qx)<U"(Qx)fla.  (B.ll) 

Thus,  we  have  shown  that  if  the  maximum  occurs  at  any  point  (?,  on  x,  =  0 ,  then 


p2  <a~2[f'\Qx)\2  +  [f\Qx)\2 

<  max  {a“2[/"(x2)]2  + 1  /'(*2)]2} , 

jr2€(0,AJ 


where  the  last  inequality  takes  (B.4)  into  account.  Taking  the  limit  as  a  ->  n/h  ,  we 
arrive  at  (3.4).  In  our  arguments  we  have  used  the  obvious  fact  that  if  the  maximum 
value  of  p1e2ax'  occurs  at  a  point  on  xx  =  0  then  the  computed  bound  for  p2e2ax' 
is  also  a  bound  for  p2 . 

It  is  perhaps  worth  mentioning  the  curious  fact  that  if  the  portion  of  the  boundary 
on  which  nonzero  data  are  given  were  not  the  flat  end  x,  =  0  but  rather  a  curve  that 
is  convex  outward  to  the  left,  e.g.,  a  semicircular  cap,  then  bounds  for  p2  could  be 
obtained  by  using  the  following  result  of  Payne  [18]: 


p  <  max 

S 


A 

K2(s) 


+f 


%S 


(B-13) 


where  f  s ,  f  Sjj  are  the  first  and  second  tangential  derivatives  of  /  on  the  convex 
curve,  respectively,  and  k(s)  is  the  curvature.  For  a  flat  end,  however,  k(s)  =  0 , 
and  the  bound  becomes  meaningless. 


Appendix  C.  Verification  of  (2.46).  In  the  motivation  for  assumption  (1.20)  given 
in  Sec.  1,  it  was  pointed  out  that  one  way  of  looking  at  this  hypothesis  was  to  think 
of  the  first  term  on  the  right  of  ( 1 .20)  as  a  measure  of  the  difference  in  p(s ) ,  p( s) 
(s  >  0)  at  different  values  of  their  arguments  and  the  second  term  as  a  measure  of 
the  difference  between  p(s)  and  p(s)  for  the  same  value  of  their  arguments.  With 
this  in  mind,  suppose  that  p  and  p  are  such  that 

|/?(s)  -  £(j)|  <  c3c2  max[p(s) ,  p(s)]  (C.l) 

for  some  constant  c3  >  0.  For  a  given  p ,  p ,  it  is  desirable  to  choose  c3  as  small 
as  possible.  (We  recall  from  Sec.  1  a  similar  observation  concerning  the  choices  of 
ma  ,  Mn,  ma,  Ma  .)  It  follows  that  if  p(s)  >  p{s) ,  then 

p(s)<P(s)<P(s)(l-e2c3)-'-  (C.2) 

and  if  p{s)  <  p(s) ,  then 

/>(s)<£U)  </>(*)(! -e2c3)-1,  (C.3) 

where  we  assume  that  1  -  e1cJ  >  0 .  Thus,  in  general 

p(s)<p(s){  l-e2c3r'  (C.4) 
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P(s)  <  p(5)(l  -e2Cj)  ‘  . 

We  look  now  at  Case  1  and  Case  2  separately. 

Case  1.  From  (C.4),  (1.16),  and  (C.5)  we  have 

p{s)  <  {A/,  +  ktsp(s)}{\  -  e2c3)_1 

<  {Ml  +  kiS[\p(s)  -  p(5)|  +  />(*)]}(1  - 

<  A/,(  1  -  e2c3)_1  +  A,1i{e2c3p(i)(l  -e2c3)“‘  +  />($)}(  1  -c2c3) 
=  A/,(l  -£2C3)_1  +  K,(1  -  62C3)~25p(5). 

Similarly,  making  use  of  (C.5),  (1.7),  and  (C.4)  we  have 

p(j)  <  A/,(l  -£2c3)~‘  +  *,(  1  - £2c3r2s/>(j) . 

It  is  also  easily  seen  from  (C.4),  (C.5)  and  (1.16),  (1.17)  that 


p{s)  >  m,(l  -e  c}), 
fi(s)  >  m,( l -e2c3). 

Now  clearly  from  (C.7),  (C.9),  and  (1.16)  it  follows  that  we  can  take 


Af,  <  Af , (I  —  <:2c3 ) _  1 ,  m,  >  mx{\ -  e2c}) , 


and  thus  we  have 


provided 


it  m.  n  m.  2  ,2  ir(l -c,)m. 

V'~  hMx~  hM ,  E  >  V'  ~  2h{  \  +  c,)A/,  ’ 


(1-ci)<2(1+c,)(1-£2c3)2. 


The  inequality  (C.12)  clearly  holds  for  sufficiently  small  e . 

Case  2.  We  now  rewrite  (C.4)  and  (C.5)  as 

tfwr'  <[p(j)(i-e2cj)r‘. 

[p{s)fl  <[p{s){l-e2c3)V' . 

Then,  from  (C.13),  (1.8),  and  (C.14)  we  obtain 

[p(s)]~l  <[M2  +  K2sp(s)]{l  -e2Ci)-' 

<  A/2(l -e2c3)_1  +K2(  1  -e2c3r2sp(s). 

The  last  inequality  in  (C.15)  follows  as  in  the  proof  of  (C.6).  Similarly, 
(1.17),  and  (C.13)  we  obtain 


(p(5)]  1  <  M2(l  -  e2Cj)  1 +*2(1  -e2c3)“25^(5). 


(C.  10) 

(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 
from  (C.14), 

(C.  16) 


and  as  before 


[^(j)]  1  >m2(l -e2c3),  [p(s)fl  >m2(l -e2c3).  (C.17) 
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Thus,  again  we  arrive  at 

*2>  (Clg) 

provided  (C.12)  holds.  We  have  thus  shown  that  if  p  and  p  satisfy  (C.  1)  and  e 
is  sufficiently  small,  then  (2.46)  will  hold.  We  remark  that,  in  the  three  examples 
discussed  in  Sec.  3,  (C.l)  is  satisfied. 
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